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Abstract. Results from a large, multi-J CO, ^'^CO, and HCN line survey of Luminous Infrared 
Galaxies (LIRGs: Lir^IO^'' Lq) in the local Universe (z^O.l), complemented by CO J=4-3 up to 
J=13-12 observations from the Herschel Space Observatory (HSO), paints a new picture for the 
average conditions of the molecular gas of the most luminous of these galaxies with turbulence 
Q ■ and/or large cosmic ray (CR) energy densities Uca rather than far-UV/optical photons from 
M , star-forming sites as the dominant heating sources. Especially in ULIRGs (Lir>10^^ L©) the 
C/3 . Photon Dominated Regions (PDRs) can encompass at most ~few% of their molecular gas mass 
while the large Ucr~10^ Ucr, Galaxy, and the strong turbulence in these merger/starbursts, can 
volumetrically heat much of their molecular gas to Tkin~(100-200) K, unhindered by the high 
y—^ • dust extinctions. Moreover the strong supersonic turbulence in ULIRGs relocates much of their 
^ ■ molecular gas at much higher average densities (^I0''cm~^) than in isolated spirals ('^(lO^- 
OO . 10'')cm~^). This renders low-J CO lines incapable of constraining the properties of the bulk 
CO ■ of the molecular gas in ULIRGs, with substantial and systematic underestimates of its mass 
possible when only such lines are used. Finally a comparative study of multi-J HCN lines and 
^ '. CO SLEDs from J=l-0 up to J=13-12 of NGC6240 and Arp 193 offers a clear example of 
■ two merger/starbursts whose similar low-J CO SLEDs, and Lir/Lco,i-o, Lhcn,i-o/Lco,i-o 
T-H ' ratios (proxies of the so-called SF efflciency and dense gas mass fraction) , yield no indications 
CN| I about their strongly diverging CO SLEDs beyond J=4-3, and ultimately the different physical 
' conditions in their molecular ISM. The much larger sensitivity of ALMA and its excellent site 
^ , in the Atacama desert now allows the observations necessary to assess the dominant energy 
sources of the molecular gas and its mass in LIRGs without depending on the low-J CO lines. 
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1. Introduction 

The CO rotational lines are the most widely used probes of the average conditions and 
mass of the molecular gas in galaxies, with a substantial body of data assembled over 
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the last two decades (e.g. Braine & Combes 1992; Aalto et al. 1995; Solomon et al. 1997, 
Downes & Solomon 1998, Mauersberger et al. 1999 Yao et al. 2003, Mao et al. 2011). 
Most of these are for J=l-0, 2-1, with only few datasets having 3=3-2 or higher-J lines. 
The lowest J=l-0 line serves as a global molecular gas mass tracer via the so-called 
Xco=M(H2)/L^o i_o factor (Dickman et al. 1986; Solomon & Barrett 1991). The low 
excitation requirements of this transition (Ei/kB~5.5K and ncrit^400cm~'^) allow even 
the coldest and lowest density gas in ordinary Giant Molecular Clouds (GMCs) to have 
a substantial, and hopefully calibratable, line luminosity contribution. 

The next two decades were then spend formulating the dependance of Xco on the 
average molecular gas conditions, and calibrating its values for various ISM environments 
(Maloney & Black 1988, Young & Scoville 1991, Wolfire et al. 1993, Bryant & Scoville 
1996. IDownes fc Solomon 19 98). Aside from a strong dependance of Xco on metallicity 
that could leave large reservoirs of molecular gas in the outer parts of spirals or metal- 
poor dwarfs untraceable by CO (Papadopoulos et al. 2002, Wolfire et al. 2010), any 
deviations of Xco from its Galactic value of ~5M0 (Kkms~^ pc^)~^ were determined 
using radiative transfer models of CO lines to constrain the average gas density n(H2), CO 
J=l-0 brightness temperature Tb,i_o, and dynamical state upon which X^o depends as 



M(H2) 3.25 Vn(H5^_i f \ 
Lco,i-o V" Tb,i-o VKkms ^ pc^ 

(Papadopoulos et al. 2012a) where a=0.55-2.4 depending on the cloud density profile, 
and Lco1-o=/av/a Tb,i-odadV is the velocity/area-integrated CO J=l-0 brightness 
temperature at the reference frame of the source. The parameter Kvir is given by 



(dV/dR) ^ ,, [C0/H2] f n(H2) V^^^ 
^^'^ - (dV/dR),,,,,, ^ ^-'^ V^A,o Wcm-^J ' ^'-'^ 

where [CO/Hal-lO-* is the CO abundance and Aco=[C0/H2]/(dV/dR) is one of the 
three parameters (the other two being n(H2) and Tkin) of one-phase Large Velocity Gradi- 
ent (LVG) radiative transfer models (with dV/dR the gas velocity gradient). The values of 
Kvir determine the average gas dynamical state, with Kvir^l-2 typical of self-gravitating 
(or nearly so) states, and KvhO^l indicating unbound gas ( [Papadopoulos et al. 2012a| . 
In principle LVG models of the lowest three CO lines can produce constrains on Xco 
to within factors of ~2 via Equation 1.1 (or via similar expressions in the literature, 
e.g. [Bryant fc Scoville 1996 1 since the excitation characteristics of J=3-2 (E3/kB=33K, 



ncrit~10^cm ■^) bracket the conditions expected for much of the gas in GMCs (Tkin^(10- 
30) K, n^ (500-10"^) cm^'^). For merger/star bursts such studies yielded alow Xco^lM© (Kkms^^ pc^)^^ 
([Solomon et al. 1997| IDownes fc Solomon 19981 lYao et al. 2003]) . which was then widely 
used in the literature for similar galaxies at low or high redshifts. 

The Xco dependance on the physical conditions of GMCs brings forth the question of 
what determines them (and thus Xco) in LIRGs. When it comes to the all-important ther- 
mal state, far-UV/optical photons from star-forming (SF) sites have long been considered 
as the main heating source of the molecular gas and dust in the Galaxy, with PDRs con- 
taining ^90% of its molecular gas in Photon Dominated Regions (PDRs) where photons 
determine chemistry and thermal balance (HoUenbach fc Tielens 1999). The molecular 
ISM and its corresponding line and dust continuum in LIRGs is then considered fully 
reducible to PDR ensembles (Wolfire et al. 1990). 
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Figure 1. The SPIRE/FTS spectra of Arpl93, NGC6240. The lines are: CO J=4-3 up to 
J=13-12 (red), the two fine structure hues of [ClfPi^^Po and ^P2-^^Pi and [Nil] (green). 



2. The line survey of LIRGs: reaching out to the warm and dense gas 

We untertook a large molecular line survey of LIRGs, drawn from the IRAS Bright 
Galaxy Survey (BGS) (Soifer et al. 1987; Sanders et al. 2003) using the James Clerk 
Maxwell Telescope (JCMT) on Mauna Kea in Hawaii, and the IRAM 30-m telescope 
on Pico Veleta in Spain. Our line database, when augmented by data existing in the 
literature, yielded the largest multi-J dataset of total GO, ^^GO and HCN line luminosities 



for galaxies in the local (z^O.l) Universe ( Papadopoulos et al. 2012b). The lines are CO 
J=1~0, 2-1, 3-2, i^CO J=l-0 and/or J=2-l, as well as HCN J=l-0, 3-2, and 4-3, while 
for a smaller subsample CO J=4-3 and J=6-5 are also available. This ground-based effort 
has been complemented by SPIRE/FTS observations of a large number of ULIRGs in 
our sample, yielding complete CO SLEDs from J=l-0 up to J=13-12, and an unparalled 
view of the densest and warmest gas in these merger/starbursts (see Figure 1). 

Even a casual inspection of the high- J CO SLEDs of the two merger/starbursts Arp 193 
and NGC 6240 (Figure 1) reveals two systems that, despite their similar low-J CO SLEDs 
and Lir/Lco,i-Oi Lhcn,i-o/Lco,i-o ratios (proxies of the so-called SF efficiency and 
dense gas mass fraction), they have strongly divergent CO SLEDs above J=4-3, indicat- 
ing different properties and/or mass for their warm and dense gas phase. This high- J CO 
SLED divergence appears even more clearly in Figure 2 where they are shown normalized 
by the corresponding IR luminosities. It is worth noting that the 00/^*^00 J=2-l line 
ratios for these two galaxies are also similar QPapadopoulos et al. 2012bl ). 

Thus the entire low-J CO, ^^CO line diagnostic that is typically used to constrain 
the average molecular gas properties (and the corresponding Xco factor) in LIRGs may 
be inadequate for the merger /starhurst systems among them. This is not unexpected 
since, unlike in isolated SF disk galaxies, the much stronger supersonic turbulence in 
the molecular gas of ULIRGs (e.g. IDownes fc Solomon 1998|) will relocate most of the 
resident molecular gas mass at average densities of n^lO^'cm"'^, and thus beyond the 
reach of the low-J CO, ^'^CO lines. The latter then can no longer yield relevant corrections 
for the Xco factor in these systems ( [Papadopoulos et al. 2012a| ). 

For Arp 193 and NGC 6240, the only ground-based observations able to break the 
aforementioned degeneracies were multi-J HCN line observations (Papadopoulos 2007). 
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Figure 2. The CO SLEDs of Arp 193 and NGC6240 normalized by their Lir (Lq/Lq units). 



These found a very low HCN(4-3) / (1-0) ratio in Arp 193, compatible with a total lack of 
gas at n^lO^cm"^, quite unlike the state of the gas in NGC6240, and what is generaly 
expected in such extreme merger /starbursts (Gao & Solomon 2004). These diferences 
are seen clearly in Figure 3 that shows the [n(H2),Tk] space compatible with the HON 
line ratios of these two systems. In the past multi-J measurements of such high-dipole 
moment molecules were impractical because of the limited sensitivities of the available 
mm/submm telescopes. This will no longer be the case with ALMA, and multi-J HCN, 
HCO+, CS observations of (U)LIRGs may be the only way of assessing the average con- 
ditions of their molecular gas (and its mass) without the aforementioned degeneracies of 



the low- J CO lines ( [Papadopoulos et al. 2012a[ ), and without depending on fully-sampled 
high-J CO SLEDs. These will remain difficult to obtain past J=6-5, even with ALMA, 
because of the large atmospheric absorption at z^^690GHz. 



3. New power sources for the molecular gas in ULIRGs 

In the high-density and metal- rich molecular ISM of ULIRGs neither far-UV/optical 
photons nor SNR-induced shocks travel far, the latter dissipating strongly in high density 
gas. The strong confinement of warm PDRs and SNR-shocked regions is evident even 
in the Galaxy where most of the molecular gas is at much lower average densities. Such 
regions involve only '-^(0. 1-1)% of the mass of atypical CMC (Papadopoulos et al. 2012b), 
and this is indeed why even the template SF clouds of Orion A, B are globally cold with 
CO(2-1)/(1-0) ratios of ^0.6-0.7, and only few isolated regions near HII regions reaching 
-0.9-1.3 (Sakamoto et al. 1994). The corresponding global CO(3-2)/(l-0) for Orion 
would then be —0.30 while for the SF "hot-spots" r32— 0.9-1.3. Thus it is indeed a surprise 
that entire galaxies can approach and even surpass these high-excitation regimes (see 
Figure 4). These excitation outliers are typically merger/starburst LIRGs, and analysis 
of their CO ratios (often supplemented by available CO J=4-3, 6-5 lines) find large 
(5^10%) fractions of dense (— (lO'^-lO^) cm~'^) and warm (—(100-200) K) molecular gas. 

For NGC 6240 and Arp 193 where HCN multi-J and fully sampled CO SLEDs from 
J=l-0 up to J=13-12 are available (Figure 1) a superposition of states from the HCN- 
constrained LVG solutions (Figure 3) with Tkin5^70K is adequate to reproduce the entire 
CO SLEDs up to the highest transition. Only for the lowest two CO lines there are 
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Figure 3. The two-dimensional probability density functions of the [n, Tkin] LVG solutions 
in steps of 0.2, as constrained by the HCN line ratios of NGC 6240 and Arp 193. Color: the 
corresponding (dV/dR) within the Kvir=0.5-2 range (Equation 1.2). 




Figure 4. The distributions of the CO (2-l)/(l-0) (3-2)/(l-0) brightness temperature ratios 
for the galaxies in the sample. The shaded area marks those for sources with CO emission sizes 
^15" (the JCMT beam). The horizontal bars indicate the ratios for the Orion SF "hot-spots". 



additional contributions from unbound low-density gas that contains little mass. For 
NGC 6240 in particular the HCN-bright phases that can reproduce its luminous high- J 
CO SLED with its large line-continum contrast (Figs 1, 2) have n(H2)^(few)x 10'* cm^"^, 
and contain most of its molecular gas mass. Far-UV photons cannot drive such high 
temperatures since for photoelectric heating it is: rpo<cAco+AH2+Aoi+Agas-dust, be- 
cause the line cooling terms (i.e. CO, H2 and 01 cooling) become large at high densities 
(AiinoOc[n(H2)]^). On other hand the strong supersonic turbulence and the high CR en- 
ergy densities in the ISM environment of (U)LIRGs like NGC 6240 can easily maintain 
large amounts of warm and dense molecular gas mass ( Papadopoulos et al. 2012a I. 
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Figure 5. The Xco factor for NGC 6240 and the conditions shown in Figure 3. 

4. A Galactic Xco in ULIRGs? 

Using Equation 1.1 with inputs from the HCN-constrained LVG solution space of 
NGC 6240 shown in Figure 3 (which also satisfies the multi-J CS and HCO+ available 
for this LIRG) we obtain the Xco factor for the prevailing conditions of its molecular gas 
(Figure 5). It is obvious that for the high-density molecular gas of this merger/starburst 
Galactic Xco values are possible. Using the SPIRE/FTS data to select the sub-regions 
of the solution space that can reproduce also the high-J CO lines selects regions with 
Tkin^70K that still include Galactic Xco values, and point towards what may be a 



systematic understimate of molecular gas mass in ULIRGs (Papadopoulos et al. 2012a). 
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